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ABSTRACT 

We present extensive new photometry in (g',i') of the large globular cluster (GC) system around 
NGC 3311, the central cD galaxy in the Hydra cluster. Our GMOS data cover a 5'.5 field of view and 
reach a limiting magnitude i 1 ~ 26.0, about 0.5 magnitude fainter than the turnover point of the GC 
luminosity function. We find that NGC 3311 has a huge population of ~ 16, 000 GCs, closely similar to 
the prototypical "high specific frequency" Virgo giant M87. The color-magnitude distribution shows 
that the metal-poor "blue" GC sequence and the metal-richer "red" sequence are both present, with 
nearly equal numbers of clusters. Bimodal fits to the color distributions confirm that the blue sequence 
shows the same trend of progressively increasing metallicity with GC mass that has previously been 
found in many other large galaxies; the correlation we find corresponds to a scaling of GC metallicity 
with mass of Z ~ M - 6 . By contrast, the red sequence shows no change of mean metallicity with 
mass, but it shows an upward extension to much higher than normal luminosity into the UCD-like 
range, strengthening the potential connections between massive GCs and UCDs. The GC luminosity 
function, which we measure down to the turnover point at Mi ~ —8.4, also has a normal form like 
those in other giant ellipticals. Within the Hydra field, another giant elliptical NGC 3309 is sitting just 
100" from the cD NGC 3311. We use our data to solve simultaneously for the spatial structure and 
total GC populations of both galaxies at once. Their specific frequencies are 5jv(AT3311) = 12.5 ± 1.5 
and Sat (A3309) = 0.6 ± 0.4. NGC 3311 is completely dominant and entirely comparable with other 
cD-type systems such as M87 in Virgo. 

Subject headings: galaxies: elliptical — galaxies: individual (NGC 3311) — galaxies: individual (NGC 
3309) 



1. INTRODUCTION 

One of the most distinctive features of the supergiant 
cD-type galaxies found at the centers of rich environ- 
ments is the very obvious presence of huge populations of 
old-halo globular clusters (GCs), often numbering above 
10,000 GCs in a singl e system an d spanning radii up- 
wards of 100 kpc (e.g. lHarrisll2001h . These "high specific 
frequency" ("high-SV') globular cluster systems (GCSs) 
provide a platform for statistical studies of halo cluster 
properties that cannot be carried out in any other envi- 
ronment. Despite their rarity, they are therefore prime 
targets for both observations and interpretive modelling 
on the formation and evolution of GCs and their host 



galaxies. 

The GCSs in cD galaxies (most of which are also 
"Brightest Cluster Galaxies" or BCGs) have been used 
to identify a new correlation between GC luminosity 
(or mass) and color (or metallicity). While the clus- 
ter metallicity distri bution function (MDF) usually has a 
bimodal shape (e.g . IZepf fc Ashmanl 119931 : lLarsen et al.1 



2001; Micsk e et all [20061 among many others), recent 
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photometric work suggests that the more metal-poor 
GCs (those in the "blue sequence") have higher heavy- 
element abundance at progressively higher luminosity. 
This mass-metallicity relation (MMR) was first discov- 
ered in eight giant BCGs by lHarris et alj ([2006, here- 
after H06) and in t he Virgo giants M87 and NGC 4649 
(|Strader et alJ l2006L The same trend can be seen for the 
GCS arou nd the giant Sa , galax y NGC 4594 (jSpitler et all 
2006) , and Mi eske et al.l (|2006l ) find evidence for a similar 
MMR in their extensive sample of Virgo Cluster Survey 
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galaxies by combining their GC data into groups by host 
galaxy luminosity. As H06 discuss, a clear trace of this 
trend ha d already been foun d for the Fornax cD NGC 
1399 by lOstrov et all (|1998f) and iDirsch et all (|2003f) . 
who noted that the bright end of the GC population 
appeared broadly unimodal rather than bimodal. As we 
now realize, this effect comes about essentially because 
the blue and red sequences merge at the top end. 

Intriguingly, the more metal-rich "red" sequence is not 
seen to exhibit a MMR of its own in any galaxy so far. 
That is, for the more metal-rich GCs, metallicity is in- 
dependent of mass. The basic interpretation proposed 
by HOG is that the blue-sequence clusters formed pri- 
marily in dwarf-sized pregalactic clouds while the first 
stage of hierarchical merging was beginning. In these 
sites, more self-enrichment could occur within the more 
massive dwarfs where more of the enriched gas from su- 
pernova ejecta could be held back. These same higher- 
mass dwarfs would also, on average, produce the most 
massive GCs CE arris fc Pud ritz 1994). By contrast, the 
more metal-rich GCs on the red sequence formed later, 
perhaps in major mergers, and at a stage when all the 
material could be held in the much deeper potential well 
of the final big galaxy. In that stage, enric hment would 
have b een more independent of GC mass. | Mieske et all 
(2006) discuss a similar view. iRothberg et all (|2007t ) use 
semi-analytic galaxy formation models to test a differ- 
ent approach, namely the contribution of stripped dwarf 
nuclei to the MMR. They find that a MMR along the 
blue sequence appears naturally with a standard model 
prescription for star formation and enrichment rates in 
semi-analytic models. An integral part of their interpre- 
tation is that many of the massive blue GCs are actually 
the dense nuclei of stripped dwarfs (dE,N or UCD ob- 
jects). 

Still more intriguingly, at least one clear anomaly 
is already known to exist: although the MMR along 
the blue sequence has now been found unequivocally 
in many galaxies, it does not seem to occur in the 
well studied Virgo giant NGC 4472, where both the 
red and blue GC sequences run vertically. This re- 
sult has been confirmed in severa l studies and in a 
variety of photometric bandpasse s (Geisler e t alJ 19961 : 
Lee et al.1 fl9p iPuzia et al.l 119991: [Rhode fc Zepl 120011 



Strader et al. 2006: lMieske et al.ll2006[) " (|Rothberg et al.1 
20071) find that increasing the efficiency of supernova 



feedback for heating in dwarfs can produce this effect, 
leaving the blue GC sequence without a MMR. But it 
is not yet clear what underlying physical or environmen- 
tal conditions would cause this to occur. In short, the 
existence of the MMR is beginning to provide rich new 
insights into the early formation era of large galaxies. 

NGC 3311, the central cD galaxy in the Hydra clus- 
ter A1060, has long been an attractive target for GCS 
studies because it is relatively nearby (d = 53 Mpc for 
a Hydra redshift cz ~ 3900 km s" 1 and Ho = 70 km 
s _1 Mpc -1 ) and has an enormously populous GCS com- 
parable with the largest known, such as M87 in Virgo 
or NGC 4874 in Coma. Its ranking among the high- 
s' at systems was established wi th the early photomet- 
ric survey bv lHarris et al.l (jl983). Fu rther ground-based 
imaging in the Washin gton system ([Seeker et al.l 119951 : 
McLaughl in et afl i~995) indicated that its GCS was mod- 
erately metal-rich with a hint of bimodality, and that its 



spatial distribution was similar to the shape of the cD en- 
velope light. Hub ble Space Telescop e imaging with the 
WFPC2 camera ([Brodie et al.l [2000) in the frequently 
used (V — I) color index was employed to argue that 
the MDF for its clusters was more like the "standard" 
pattern in giant ellipticals - that is, a probable bimodal 
structure with metal-poor and metal-rich sequences at 
their usual mean metallicities. T he most recently pub- 
lished study ([Hempel et al.ll2005l ) combined the WFPC2 
(V, I) data with NICMOS imaging in F160W (approx- 
imately the H— band) to gain a first rough estimate of 
the metallicity and age distributions through models in 
the (V — I, V — H) two-color diagram. 

None of these earlier studies have given a satisfac- 
tory picture of the system. The HST -based studies 
([Brodie et al.l 120001 : iHempel et all l2005f ) are hampered 
by the small fields of view of WFPC2 and NICMOS, 
yielding very incomplete spatial coverage of this big cD 
galaxy, and thus relatively small G C samples to work 
with. By contrast, the initial study of lHarris et al.l (j!983l ) 
in the pre-CCD photographic era covered a wide field 
but penetrated to only quite shallow photometric lim- 
its by today's standards. Relatively shallow limits and 
photometric calibrati on problems also a ffected the first 
CCD-based imaging ([Seeker et al.lll995ft . To study this 
unusually attractive GCS further, we need to draw from 
a new photometric study that combines depth with wide 
field. 

2. OBSERVATIONS AND DATA REDUCTION 

As part of a new program (Gemini GS-2006A-Q-24) 
intended to obtain deep multicolor imaging of globular 
cluster systems around cD galaxies, we obtained (g', i') 
images of NGC 3311 and the central Hydra cluster field 
using the GMOS imager on Gemini South, which has a 
5.5' x 5.5' field of view (FOV) and a scale of O.U6"/pix 
after 2x2 binning. The raw observations were taken on 
the nights of 2006 February 8 and March 23, under dark 
skies and photometric conditions, with an average seeing 
of 0.5". In each waveband, 13 individual exposures were 
taken, which we then reregistered and median-combined 
to create the final, deep images in (g',i'). The total 
integration time in each filter was 3900s. 

The preprocessing of the raw exposures was completed 
with the GEMINI package in IRAF 1 . In the i' frames, 
a small amount of night-sky fringing was present, which 
we successfully removed with a calibration fringe frame 
from the Gemini archives. After subtracting the fringe 
frame and then median combining all 13 i'-band frames, 
the resulting fringe signal contributes only 0.3 percent to 
the background fluctuations. 

For the photometric ca libration, we u sed exposures of 
Landolt standard stars (lLandoltl Il992h taken adjacent 
to our program-field images during photometric condi- 
tions. We transformed the Landolt standard-star indices 
in UBVRI to ( g',i') using the conversion equations of 
iFukugita et~aTI (l996h 2 . Since our standard star expo- 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 A complete set of the lLandoltl 1119920 standard stars trans- 
formed from the Johnson-Cousins system into the Sloan filter 
set has been compiled by the authors and is available online at 
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Fig. 1. — Photometric completeness of detection in g' and i' . 
Here the completeness fraction /, as determined from the artificial- 
star tests described in the text, is plotted as a function of magni- 
tude, where f(g') is plotted in the solid dots as a function of g', 
and /(«') is plotted in the open symbols as a function of i' . The 
interpolation curves drawn through each set of measurements are 
the two-parameter functions listed in the text. 



TABLE 1 
Basic Parameters for 
NGC 3311 



Parameter 



Value 



Type 
a (J2000) 
5 (J2000) 
v r (helio) 

Ay 
(m - M) 
E(g' - i<) 

v i 



E2/cD 
10 ,i 36 m 42?8 
-27?31'42" 
3593 km s" 1 
0.26 
33.62 
0.158 
11.04 
-22.8 



sures covered only a small range in airmass, we used the 
NGC 3311 images taken within just the same airmass 
range, and taken during the same night within one hour 
of the standard-star exposures. This calibration was suit- 
able only for establishing a mean zero-point value for our 
NGC 3311 frames, and not for a more complete evalua- 
tion of the coefficients in airmass or color. We therefore 
regard the calibration as only a preliminary one, to be 
refined when a larger network of standards can be ob- 
tained. 

The final measurement of our combined, deep (g 1 
images was carried out with the latest im plementation 
of th e standalone daophot code, daophot4 (jHarris et al.l 
|2007| L We used an extensive series of artificial-star tests 
with the addstar component of daophot to evaluate the 
internal photometric detection completeness as a func- 
tion of magnitude: stars were added to each the g' and 
i' images in groups of 500 over a wide range of magni- 
tudes, the images were then remeasured in the same way 
as the original frames, and the fraction / of stars recov- 
ered in each magnitude bin was then determined. The 
results are shown in Figure [TJ The limits of our data, 
defined as the magnitudes at which / drops to 0.5, are 

http: / / www.elizabethwehner.com / astro / sloan.html 



g'(lim) = 26.67 and i'(lim) = 26.30. The trend of / 
with magnitud e is well matched by a Pritchet interpola- 
tion function ([Fleming et al.|[T995l ). 
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g(m - mum) 
yj\ + a 2 (m - m Um y 



(1) 



which has two free parameters: the limiting magni- 
tude mum, and a parameter a giving the steepness of 
the dropoff. For these images we find a g i = 0.8 and 
= 1.4. As will be seen below, the limiting magnitudes 
g'(lim) = 26.7 and i'(lim) = 26.3 are deep enough to be 
comparable with the GC luminosity function turnover 
point of NGC 3311. Within the bright bulge of NGC 
3311, the photometry depth diminishes due to the in- 
creased local sky noise at small radii. However, for radii 
larger than 0.'6 from NGC 3311, which include the vast 
majority of our data, we find no change in the complete- 
ness limits g' , i'(lim). 

The random uncertainties of the photometry, as deter- 
mined from daophot /allstar and from the same artificial- 
star tests, are well approximated by the interpolation 
curves 

a(g') = 0.01 + 0.050 exp(</ - 26.0), (2) 
o-(i') = 0.01 + 0.050 exp(i' - 25.0). (3) 

They show little variation from place to place within the 
GMOS field, indicating that crowding is unimportant, 
and that background light from NGC 3311 and NGC 
3309 has little influence except near their very centers. 
On the basis of these error curves, we would expect the 
GC sequences to start being noticeably 'broadened' in 
color, purely because of random uncertainties in the pho- 
tometry, for levels i' > 24. Lastly, the artificial-star tests 
show that any systematic error (net bias) in the photome- 
try is less than 0.02 mag in each filter for i' < 25, g' < 26. 

Throughout the following discussion we adopt a dis- 
tance modulus (m - M) = 33.62 for NGC 3311 and 
foreground redde n ing E (q' — i') = 0.158, as we did in 
iWehner fc Hirrii (12001 . 

3. THE COLOR-MAGNITUDE DISTRIBUTION 

After carrying out the photometry in g' and i' 
with daophot we used t he additional code SExtractor 
(|Bertin fc Arnoutill999l ) to help reject obvious nonstel- 
lar objects from the detection lists (at 0.5" seeing and 
at the distance of the Hydra cluster, globular clusters 
appear indistinguishable from the starlike PSF. Under 
HST-type resolution of 0.1", the profiles of the very 
biggest GCs can, however, be resolved; see Wehner & 
Harris 2007). Detections from the two wavebands were 
then correlated, and only those objects found in both 
were kept for the final analysis. 

The final color-magnitude measurements for 8108 star- 
like objects across the entire GMOS field are shown in 
Figure O The globular cluster population is readily 
visible as the dominant concentration of points in the 
color range 0.6 < (</ — i') < 1.4. Fainter than i' ~ 25 
the field contamination and photometric scatter increase 
markedly, but brighter than this, the GC population 
is well defined. Evidence of the standard "blue" and 
"red" sequences can already be seen and will be dis- 
cussed more quantitatively below through objective bi- 
modal fits. On the red-sequence side (approximately in 



4 



Wehner et al. 



the range 1.0 < (<?' — i 1 ) < 1.3) a distinctive bright ex- 
tension to high luminosity can be seen that, intriguingly, 
does not exist for the blue GC sequence. This bright-end 
feature extends up into the GC mass range M > 10 7 Af Q 
that overlaps with the regime of the "Ultra-Compact 
Dwarfs" (UCDs). We have discussed thes e objects in 
our previous paper (|Wehner fe Harris! l200l . 

To date, not much other photometry of globular cluster 
systems in the (g' , i') indices has been published that we 
can make comparisons with. iForbes et al.l (|2004l ) present 
(g 1 — i 1 ) data taken with the GMOS camera on Gemini 
North for the Virgo giant NGC 4649. The bimodal fit 
they obtain to its GC color distribution has mean val- 
ues (</ — i')o for the blue and red sequences which are 
~ 0.1 mag redder than our calibrated and dereddened 
colors for NGC 3311. We would expect a priori, how- 
ever, that the colors of the two sequences w ould be quite 
similar from one galaxy to the next (e.g. lHarris et al.l 
120061 IStrader et aLll2004h with only a modest trend due 
to the parent galaxy luminosity. For the present, we do 
not regard this difference as significant, because (as em- 
phasized above) we regard our zeropoint calibration as 
a preliminary one, which still needs a separate series of 
observations especially to establish the photometric color 
terms. Notably, Forbes et al. also regard their own zero- 
point calibration as preliminary. What is perhaps more 
important at this stage, and equally worth noting, is that 
the mean color difference A(g' — i') ~ 0.30±0.08 between 
the red and blue sequences that we find for NGC 3311 
(sec the discussion below) is identical with the separation 
btween those two s equences measured for NGC 4649 by 
IForbes et all (120041 ). 

4. THE GLOBAL FEATURES OF THE GLOBULAR 
CLUSTER SYSTEM 

In sections 4.1 - 4.2 we use our GMOS data to discuss 
the spatial distribution of the GCs around both NGC 
3311 and 3309, including the radial and azimuthal com- 
ponents. A continuing theme in this analysis is the rel- 
ative contribution of the two galaxies to the GC popu- 
lation we see on our images. In section 4.3, we discuss 
the color (metallicity) distribution of the GCs by radius 
around NGC 3311, and in sections 4.4 - 4.5 we discuss 
the GC luminosity function and specific frequencies. 

4.1. Measuring the Radial Profiles: The NGC 3309 
Problem 

The core of the Hydra cluster contains several other 
large galaxies in addition to the central supergiant that 
was our primary target. Of these others, the closest in 
projection and by far the most luminous is the giant el- 
liptical NGC 3309, just 100" away from the center of 
NGC 3311. Its integrated magnitude (Vr = 11.94; see 
below) places it 0.9 magnitude fainter than NGC 3311 
(Vr = 11.04), though the true total luminosity of NGC 
3311 is estimated to be about 1.3 mag brighter than that 
of NGC 3309 because of its extended cD envelope. Nev- 
ertheless, we might expect that NGC 3309 would con- 
tribute significant numbers of GCs of its own to the popu- 
lation of objects in the field. By using starcounts around 
each of the two galaxies on the sides opposite to th e other 
galaxy. lHarris et alJ (|1983h and lMcLaughlin et all (|1995l ) 
both argued that NGC 3309 was likely to be contribut- 
ing no more than a few percent of the GCs in the Hydra 



core field, with NGC 3311 acting as the dominant con- 
tributor. Even though NGC 3311 is clearly a high-S^r 
elliptical like many cDs, such claims might mean that 
NGC 3309 would have unusually low SV for a giant el- 
liptical in a rich cluster. A better and more quantitative 
estimate of the relative populations, as well as the spatial 
structures of the two GCSs, is possible from our deeper 
data. 

To quantify our analysis, let us assume that the pro- 
jected number density of GCs on the sky from each of 
the two galaxies can be modelled by a circularly sym- 
metric function er(r) giving the number of GCs per unit 
projected area (we justify the assumption of symmetry 
below). Then, the total number density of counted ob- 
jects at a point (x, y) anywhere in the field is given by 



a{x,y) = a b + eri(ri) + cr 2 (r 2 ) 



(4) 



where r\ is the distance of the given point (x, y) from the 
center of NGC 3311, r 2 is its distance from the center of 
NGC 3309, and tj\, 0-2 are the two radial functions spec- 
ifying the GCS around each galaxy. Finally, a b is the 
uniform background density of whatever contaminating 
objects fall within our selected sample. Our goal is to 
deduce, from a fit to the observed GC distribution in the 
field, the best-estimate parameters for o\ and ct 2 simul- 
taneously. 

To perform the fit, we adopt two possible models for 
the distribution: 
(a) A standard Hubble profile, 



<r(r) = 6(1 + -) " 



(5) 



which is a power law with an added core profile of core 
radius c. 

(b) A generalized Sersic profile, 



a(r) = a e exp -b^-) 1 ^ - l) 



(6) 



where b n ~ 1.99n-0.3271 (|Ferrarese et al.ll2006D to guar- 
antee that r e will be the half-light radius, and for typical 
E galaxies, the exponent n is found observationally to 
fall in the range ~ 1 — 4. The case n = 4 recovers the 
standard de Vaucouleurs law. 

To carry out the actual numerical fit, we divide up 
the observed area into k x k boxes, each of area A; we 
use the model parameters to predict how many objects 
Nij = A-a(xi,yi) should be found in each box at position 
(xi, yj), i, j = 1, ...fc; and then calculate the residuals (ob- 
served — predicted) in each box. A standard goodness 
of fit based on Poisson statistics can then be formed, 



x 2 = 



E 



(Njjjobs) -Njjjpred))' 
Nijipred) 



(7) 



The outline of the model is shown in Figure [3] for the 
14 x 14 grid that we used for the final reductions. An ad- 
vantage of this procedure is that it uses the information 
everywhere in the frame, not just the counts within some 
arbitrary radius centered on NGC 3311. In practice, we 
ignored any boxes located at the centers of the two galax- 
ies and at one very bright foreground star, as well as the 
boxes under the shadow of the camera guide probe (lower 
right part of Fig. [3]) and at the vignetted corners of the 
field. To reduce field contamination from the outset (see 
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Fig. 2. — The color-magnitude diagram for all measured objects 
in the CMOS field centered on NGC 3311. 
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Fig. 3. — Grid layout showing the determination of the spatial 
structure of the GCS. Here North is at left and East at top, fol- 
lowing the instrumental orientation of the GMOS camera during 
the observations. The center of NGC 3311 is near the center of 
the field, while the center of NGC 3309 is at lower left; the guide 
probe shadow is at lower right. The dots are the measured objects 
in the restricted range (i 1 < 25.0, 0.6 < (g' — i') < 1.4 as described 
in the text. Concentric circles are drawn around each of the two 
major galaxies, while the grid lines show the 14 X 14 matrix of 
boxes used to solve simultaneously for the radial distributions of 
the GCSs around both galaxies. 

Fig. [2]), we used the sample of 4393 objects within the 
restricted range i' < 25.4, 0.6 < (g' — i') < 1.4, expecting 
that this sample would be dominated by GCs. 

When we fit Eq.[4]to the data, in principle we have a to- 
tal of six free parameters. For the Hubble profile model, 
these six are aj>, a\, c%, a^, C2, and the relative amplitude 
&i/&2- For the Sersic model, they are <7&,r e i,r e 2,ni,n2 
and <7 e i/<7 e 2- Once a particular ratio &1/62 is chosen, 
the actual values of both b\ and hi individually are au- 
tomatically set by the condition that the totai predicted 



population of objects in the field must precisely equal the 
observed total. The same applies to a e \ and a e i for the 
Sersic model. 

In practice, the background density ah is extremely 
hard to determine by the fitting procedure alone, be- 
cause the NGC 3311 GCS is so extensive that it domi- 
nates the starcounts even near the edges of our GMOS 
field. An external constraint on a\> is therefore desirable. 
Lacking any remote "control field" to rely on, we use the 
distribution of stars in the CMD itself (Fig. [5]) to esti- 
mate the background. The background contamination 
will be due to a combination of foreground Milky Way 
stars and faint, small remote galaxies. We can estimate 
the number and color distribution of possible field stars 
us ing the gala c tic st ellar population synthesis models 3 
of iRobin et aD (|2003h . Using a 5.5' x 5.5' field centered 
on NGC 3311, the models find only 13 stars, suggest- 
ing contamination from galactic foreground stars is small 
at these latitudes. Using a larger field of view, such as 
a square degree, results in over 2 thousand stars with 
a color distribution leaning blueward of the blue peak 
of the GCS distribution. Our main source of contami- 
nation, as indicated by the initial FOV search, is most 
likely distant background galaxies. Without a control 
field, we cannot specifically quantify these background 
galaxies. If we assume that their distribution of colors is 
roughly uniform over the interval 0.0 < (g' — i') < 2.0, 
we can use the numbers of objects with colors in the 
ranges 0.0 < (g' - i') < 0.6 and 1.4 < (g' - i') < 2.0 
to predict how many background objects are in our se- 
lected midrange of 0.6 < (g 1 ~ i') < 1.4. For i' < 25.4 
we count 860 objects total that are bluer or redder than 
the GC color range. Normalizing to the same color inter- 
val, we then estimate that 570 of the 4393 objects within 
0.6 < (g' — i') < 1.4 (or about 13%) are background 
contamination, equivalent to o\, — 18 arcmin~ 2 . In the 
numerical fitting procedure we adopted this background 
throughout and treated only the 5 remaining parameters 
in each model as free. 

For the Hubble profile model, we obtain a best-fit set 
of parameters a x = -2.04 ± 0.07 and a = 94" ± 8" 
for NGC 3311; and a 2 = -2.8(+0.7, -1.3) and c 2 = 
55"(+30", -24") for NGC 3309. The ratio (&2/&1) is 
0.093±0.057. For this set of parameters the minimum re- 
duced Xv value was 1.44. We experimented with different 
grid sizes and found that the results were not sensitive 
within broad limits to the fineness k 2 of the chosen grid. 
Given the simple nature of the model, and other factors 
not accounted for such as the necessarily patchy sam- 
pling across the grid, the internally uncertain adopted 
level of the background, and possible dumpiness (clus- 
tering) in the background population, the fit is successful 
at accomplishing our main goal to evaluate the relative 
contributions of the two giant galaxies. 

A very slightly better overall fit (minimum \v = 1.41) 
is obtained with the generalized Sersic model. For the 
NGC 3311 GCS, we find m = 1.24 ± 0.06 and r el = 
177" ± 9", while for NGC 3309 n 2 = 1.1(+1.2, -1) and 
r E 2 = 93" (+78", -45"). The ratio of the two amplitudes 
is (a e2 /a el ) = 0.080 ± 0.046. 

For comparison, the scale sizes of the underlying 
halo light of each galaxy are significantly smaller: 

3 http:/ /bison. obs-besancon.fr/modele/ 
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Fig. 4. — Radial distribution of the GCs in the field relative to 
the center of NGC 33f 1. The number of objects within each of 
the grid boxes shown in the previous figure is plotted against the 
radial distance of the box center from NGC 3311. The best-fitting 
Hubble profile (see text) is shown as the dashed line, while the 
best-fitting generalized Sersic model is the solid line. 

IVasterberg etHI (fl99iT) find r e ~ 21% ~ 5.6 kpc for 
NGC 3309 and 94'.'2 ~ 25 kpc for NGC 3311 from an 
integrated-light study in the Gunn r filter. These effec- 
tive radii are not strictly comparable to our solutions for 
the GCSs, because they assume a classic de Vaucouleurs 
profile (n — 4 in the Sersic model); but they confirm that 
the GCS is spatially mo re extended, as has usually bee n 
found for giant galaxies (iMcLaughlinll 19991 : lHarrisfeoOlf ) . 

One very evident result of this numerical exercise is 
that NGC 3309 contributes only a minor part of the GC 
population in our field (2 to 3 percent; see section 4.5 
below). As a consequence, its GCS shape parameters are 
quite weakly constrained. Assuming that the detected 
objects in the field belong entirely to NGC 3311 and 
ignoring NGC 3309 completely (that is, setting er 2 = 
in the numerical fitting) yields a minimum \ 2 negligibly 
worse than the two-galaxy fit. We quantify the relative 
importance of the two galaxies a bit further below. 

A rough comparison of the two model fits is shown in 
Figure 21 where we simply plot the number density of 
objects in each grid box as a function of distance from 
the center of NGC 3311. The best-fit Hubble and Ser- 
sic models for NGC 3311 alone are shown overplotted, 
verifying the first-order description of the entire popu- 
lation in the field as belonging to the central cD. The 
same grid plot, but this time centered on NGC 3309, 
gives a radically different impression, similar to picking 
just a random point somewhere away from NGC 3311. 
That is, no significant concentration towards NGC 3309 
is evident. 

A significant improvement in this estimate of the two 
GCS spatial structures will be possible if a better ex- 
ternal measurement of the background can be obtained. 
The NGC 3311 system is so extensive, however, that pho- 
tometry at projected radii r > 10' (150 kpc) will be nec- 
essary for a control-field location. 

4.2. Azimuthal Distributions 

McLaug hlin et all (|1995l ) also attempted to measure 
the ellipticity of the NGC 3311 GCS and to compare it 
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Fig. 5. — Azimuthal distribution of the GCS around NGC 3311. 
The numbers of objects within 30-degree sectors, and within radii 
0^26 < r < 2'.50 of NGC 3311, are shown as the data points with 
error bars. The total sample includes 3544 objects divided into 30- 
degree sectors; the two open circles denote the sectors containing 
the center of NGC 3309 and the shadow of the GMOS guide probe. 
These are not used in the numerical fit. The sinusoid-like solid 
curve (see the text for the actual model curve) is the best-fitting 
ellipsoidal model, with e = 0.07 and major axis 9 p = 90°. For 
comparison, the solution obtained by McLaughlin et al. (1995) from 
a smaller photometric sample and a brighter limit is shown by the 
crosses with error bars and the dashed curve. 

with the shape of its cD envelope. We can, in princi- 
ple, make a better test with our much larger measured 
population. In this case, for our test sample we take all 
3544 measured objects within radii 0.'26 < r < 2'.50 of 
NGC 3311 (the outer circle is the largest complete one 
fitting within the GMOS frame boundary) and within 
the photometric limits i' < 25.4,0.6 < (g' - i') < 1.4. 
This sample is subdivided into 30° sectors and plotted 
against azimuthal angle 8 of the sector, as shown in Fig- 
ured] The two sectors containing the center of NGC 3309 
and the guide probe shadow are shown as open symbols 
in Fig. [5] and are ignored. 

It is evident that there is no strong trace of ellip- 
ticity in the residuals about the mean level of ~ 275 
objects per sector after background subtraction. Al- 
though there are sector-to-sector differences larger than 
the internal statistical uncertainties, they cannot easily 
be made to match an approximate double-cosine curve 
that would be the mark of a modes tly elliptical distribu- 
tion (|McLaughlin et al.|[l994l . ll995h . The integrated halo 
light of NGC 3311 has a small but defin ite ellipticity of 
e = 0.2 and major axis at ^ 30°/210° (jDisnev fc Walll 
119771: IVasterberg et al.lll991h . If the formation of the 
GCS is closely associated with the galaxy halo, we should 
expect no large differences in their intrinsic shapes and 
orientations. Contrarily, a dramatic difference in the 
GCS shape or orientation relative to the halo light would 
be a signal of a more fundamental difference in the GCS 
formation era or its internal dynamical evolution. 

We have carried out a solution for e and 9 V with 
our GCS da ta, using the formalism of iMcLaughlin et all 
(|1994fl995l ). They show that for an intrinsically ellipti- 
cal distribution of points sampled in circular annuli, the 
observed number density is 

<r(r,0) = a b + a r- a [cos 2 {8-8 p )+(l-e)- 2 sm 2 (6-e p )]~ c 

where the radial falloff is modelled as a simple power law 
of exponent a. Here, 6* p ±180° denotes the position angle 
of the major axis at which o~(r, 8) reaches its maximum 
value. For NGC 3311, we have a ~ 1.2 (for the Hubble 
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profile fit described above, we found a = 2.0, but if we 
remove the core-radius parameter c, we require a flatter 
overall power law to match the whole profile adequately) . 
We then vary the assumed e and 8 p till the residuals 
about this model are minimized. 

The best-fit solution we find is for e = 
0.07(+0.11, -0.07) and 6 p = 90° ± 50°, shown as 
the solid curve in Fig. [5] The corresponding major axis 
line is thus 90°/270° ± 50°. (In Fig. [31 the major axis 
would run vertically.) The best-fit solution is weakly 
determined with a shallow minimum at \v — 1.95, 
and indeed solutions with e ~ are scarcely worse 
than this one. Our deduced ellipticity is lower than 
for the integrated light, and our solution for Q p falls 
~ 60° eastward of the integrated halo light, though 
its uncertainty is large enough that there is no formal 
disagreement. In short, we appear to be confronted with 
no difficulties by assuming first-order circular symmetry. 

For comparison , we show one of the solutions by 
McLaughl in et al.1 (|1995h as the dashed line in Fig. EJ 
this clearly does not match our solution very well. 
iMcLaughlin et all (|1995l ) performed solutions for two dif- 
ferent magnitude ranges (Ti < 23, T% < 24) and for the 
more slightly more restricted radial range 0.'89 < r < 
2122. For 7\ < 23, they found e = 0.32 ±0.05, p = 145°, 
while for Ti < 24 they found e ~ 0.3, 6 p = 190° (this 
latter one is shown in our Fig. [5|) . They comment that 
"the uncertainty in 9 p is at least ±40°" . 

Both of the McLaughlin et al. solutions differ from 
ours; we find a distinctly smaller ellipticity and a posi- 
tion angle between their two solutions (although, again, 
it is not clear that any genuine formal disagreement ex- 
ists). To attempt to trace down causes for the evident 
differences, we rebinned our data to magnitude ranges 
corresponding more closely to theirs, and also restricted 
to exactly the same radial range as they used. For these 
numerical trials we also included objects at all colors, 
as they did. For our objects brighter than i' = 23, we 
find e = 0.34 ± 0.15 and 6 p ~ 110°/290°. For a slightly 
deeper cut down to i' = 24, we find e = 0.10 ± 0.07 and 
9 p ~ 75°/255°. Interestingly, the best-fit position angles 
of our two subsamples bracket our best solution from the 
large, 'best' dataset [i' < 25.4 and color-selected), but re- 
stricting the sample much more severely to the brightest 
objects tends to lead to a larger e, as found by McLaugh- 
lin et al. Without being able to offer definitive reasons for 
the disagreement between these two studies, we suggest 
that the earlier GCS result was likely to have been com- 
promised by the combination of small-number statistics, 
and especially by much larger relative field contamina- 
tion (since the dominant source of field contamination 
is faint, small background galaxies which have clumpy 
space distributions, they are quite able to bias our so- 
lutions for the spatial distribution if they contribute a 
high fraction of the total sample). In the present data, 
we have been able to reject non-GC objects more com- 
pletely by color and morphology than in the older study. 

The most intriguing feature of the azimuthal distri- 
bution is the residual difference between the GC system 
(e ~ 0.1, 9 P ~ 90°) and the halo light (e ~ 0.2, 6 P ~ 30°). 
Unfortunately, we are not yet in a position to place heavy 
significance on this, simply because the solution for the 
best-fit position angle is extremely uncertain if the el- 
lipticity is very small. If the difference is verified with 
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Fig. 6. — Upper panel: Color index (g' — i')g of globular clusters 
in our CMOS field, plotted as a function of radial distance from 
NGC 3311. All objects brighter than i = 24.0 are included here. 
The overlaid diamonds indicate the RMIX fits for the blue and 
red peaks of the distribution. Lower panel: Cumulative radial dis- 
tributions of the metal-poor clusters (dashed line) and metal-rich 
clusters (solid line). The lower-metallicity clusters are slightly more 
extended spatially than the metal-richer ones, but not significantly 
so. 

deeper and cleaner samples, it will, perhaps, be a tanta- 
lizing hint that the GC accretion or destruction histories 
differed in second-order ways from the field stars. 

4.3. Radial Color Distribution 

With the knowledge that the vast majority of GCs in 
this field belong to NGC 3311, we can also investigate 
their color distribution as a function of galactocentric 
distance. In many GCSs within large ellipticals, a color 
gradient can be found in the sense that the mean color of 
the GCs (averaging together both the red and blue ones) 
conventionally decreases with increasing R ac . However, 
beginn ing especially with the key work of iGeisler et al.l 
(1996) on NGC 4472, it became clear that (in Geisler et 
al.'s concise wording) "most of this gradient appears to 
be due to the varying radial concentrations of the two 
populations". That is, the metal-poor, blue GC popula- 
tion is less centrally concentrated than the metal-richer 
population and thus contributes a higher proportion of 
clusters at larger radius, generating an overall declining 
metallicity gradient in the GC population as a whole. 
Notably, each of the two metallicity components shows 
little or no color (metallicity) gradient of its own. Very 
much the same conclusion had already become clear long 
ago (albeit with a much smaller sample of GCs than can 
be found in giant ellipticals) for the Milky Way GCS 
(lZinnlll985ft . 

In Figure[6]we show the color indices for the objects in 
the NGC 3311 field (in order to minimize field contamina- 
tion, we restrict the sample to the 1360 objects brighter 
than i' = 24.0). A broad range in color is present at all 
radii, but no strong radial gradient in color is evident for 
either the red or blue subpopulations, a similar result to 
other galaxies. In order to determine whether there was 
a trend in color versus radius, we binned the clusters in 
radius and performed a fit using RMIX (a fitting code 
discussed in detail in section 6). The locations of the red 
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and blue peaks are shown as diamonds in Figure [6j A 
comparison between the cumulative distributions of the 
red and blue clusters (shown in the lower panel of the Fig- 
ure) indicates that any difference has weak significance, 
at least within the regions our data are sampling. Un- 
fortunately, the field size we sample reaches a maximum 
R gc ~ 50 kpc, which is not large enough to test either 
the true spatial extent of the system or what happens to 
the color distribution near its outskirts (recall, for exam- 
ple, that the M87 GCS e xtends detectably ou tward to at 
least R gc ~ 100 kpc; see lTamura et al.l ((2006(1 ). 

4.4. Luminosity Function 

Deep imaging from galaxies of all types now shows 
that the luminosity distribution of their GCs follows a 
roughly universal shape that has most often been approx- 
imated by a Gaussian distribution i n number of GCs per 
unit ma gnitude (e.g. jHarrid 1200 It iKundu fe Whitmord 
I2001alr3 : lLarsen et all 120011 : Uordan et al.l 120071 among 
many others). For giant galaxies, the turnover or peak 
freq uency of the GCLF is consistently near M v = — 7.4± 
0.2 (|Harrisl[2001f ) while the intrinsic dispersion is in the 
range a = 1.3 — 1.5 (see particularly Uordan et al.|[2007l 
for an extensive analysis of the correlation of the turnover 
and dispersion with host galaxy luminosity, although it 
is worth noting that the only cD-type galaxy included in 
their data is M87). 

Ideally, a determination of the GCLF for our target 
NGC 3311 would best be done by taking the LF of an 
adjacent "control field", observed under identical condi- 
tions, and subtracting it from the LF of the objects in 
our 5.'5— wide central field. Lacking this, a second-best 
alternative is to use the outer parts of our GMOS field as 
a "background". In fact, much of the outer zone is still 
populated by GCs, and so by defining it as background 
we will partially subtract some of the GC population 
as well as the true background. However, as long as no 
large changes in the GCLF with radius exist, the residual 
GCLF we obtain this way will be systematically correct 
even if a bit oversubtracted. 

An inevitable limitation of using the outer parts of the 
image field to define the background level is that the 
background estimate <7b is much more uncertain numeri- 
cally than if we had a large, remote control field to work 
with. We divide our field into two regions: an inner an- 
nulus of inner radius ro, outer radius r\, and area Ai n \ 
and an outer zone of inner radius T\ and outer boundary 
set by the edge of the frame, and area A out . We set r 
to be the innermost radius at which the photometry is 
still highly complete to the magnitude limit of interest 
(i 1 ~ 26). The only significant free parameter we have to 
define these two zones is the mid-radius r% . The residual 
GC population within Ai n after background subtraction 
will clearly be 

N r = N in - (A m /A out )N out (9) 

where, by hypothesis, the background is constant 
across the field and subtracts cleanly out after the area 
normalization. We would like to choose the midpoint T\ 
in such a way as to minimize the relative uncertainty in 
a(N r )/N r . If n is too large, then N r will be large, but 
the area normalization factor (Ai n /A out ) will blow up 
and the background correction becomes extremely uncer- 
tain. But if n is too small, then so is Ni n , and again the 



relative uncertainty in the residual population becomes 
large. 

To fix r*i , we adopted the radial profile parameters de- 
scribed above for the GC system (the best-fit Hubble pro- 
file plus the background level crj = 18 arcmin~ 2 ). The 
inner radius ro is set at 0.'5. We then calculated the re- 
sulting relative uncertainty a{N r )/N r as a function of ri, 
where we let r\ vary from just slightly larger than ro out 
to almost the borders of the field. The optimum n turns 
out to be roughly halfway from the center of NGC 3311 
to the edge of the field, given the geometry of this system, 
but the relative uncertainty goes through a very shallow 
minimum near this point and the precise choice is not 
critical. For the following discussion, we use r\ = 2.'5. 
The "background" region is defined to be everything be- 
yond that radius, though not including the vignetted re- 
gions in the outskirts. In each zone we count all objects 
within the broad color range 0.2 < (g' — i') < 1.4, and 
also correct the raw numbers of objects in each bin for 
photometric incompleteness / = /•£'■ f g > with the com- 
pleteness function parameters derived above. 

The results are in Figure [3 shown in Ai' = 0.2— mag 
bins. For i' > 26.5, the net completeness / falls be- 
low 0.5 and the LF becomes considerably more uncertain 
(note the rapid growth of the internal error bars). This 
level is also very near the predicted magnitude of the 
GCLF turnover point for a normal giant or supergiant 
elliptical galaxy. For M v ~ -7.4 and (V - 7) = 1.0, 
we expect M° ~ -8.4 or 1° = 25.22 adding the dis- 
tance modulus (m — M)q — 33.62. To convert I to i' , 
we use the mean difference (i' — I) = 0.28 ± 0.02 deter- 
mined directl y within our field from 12 of the bright- 
est GCs fsee I Wehner fc Harris! 12007( 1 . Finally adding 
the foreground extinction A\ ~ A41 = 0.15, we have 
i' (turnover) ~ 25.65. A Gaussian curve with that peak 
magnitude and with a — 1.5 mag is shown superimposed 
on the observed GCLF in Fig. [7) Encouragingly, the raw 
LF of all objects within the inner zone (shown as the 
open circles) has the same peak to within ~ 0.1 mag 
even without background subtraction. 

The Gaussian curve shown in Fig. [7] is not intended to 
represent a rigorous derivation of the observed turnover 
point or dispersion; we use it only to show that the stan- 
dard values provide an entirely plausible description of 
the data for NGC 3311. That is, the GCLF in this cD 
galaxy gives no indication of being unusual. Our pho- 
tometry does not reach far enough past the turnover to 
permit any kind of useful test of more sophisticated nu- 
merical models, such as the asymmetric S c hecht er-likc 
function explored in depth bv I Jordan et al.l ((20071 ) from 
the Virgo cluster survey. However, in a broader sense our 
data clearly show that imaging with ground-based 8m- 
class telescopes is capable of penetrating deep into the 
globular cluster populations of galaxies in the 50-Mpc 
distance range and even beyond. 

Recent interest has followed the possibility that the 
GCLFs may differ between the redder and bluer metal- 
licity groups of clusters. If the underlying mass distribu- 
tions of both types of clusters are similar, then the GCLF 
turnover magnitude should differ between photometric 
bandpasses by typically ~ 0.2 m ag purely because o f 
differential bolometric correction ((Ashman et al.lfl995jl . 
Larger differences between the red and blue GCLFs (or 
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Fig. 7. — Luminosity function for the globular clusters in NGC 
3311. The open circles with errorbars show the numbers of objects 
within the radial range 0'5 < r < 2.'5. The connected solid dots 
with errorbars are the totals after approximate background sub- 
traction, where the CMOS image region outside the circle r > 2'.5 
is used to define the background. The red solid line is a Gaussian 
LF with peak at i' = 25.65 (corresponding to My — —7.4) and 
with standard deviation <r,/ = 1.5. 




Fig. 8. — Luminosity function for the globular clusters in NGC 
3311, plotted separately for the blue and red GC sequences. The 
open circles with errorbars show the redder clusters lying in the 
color range 0.8 < (</ — i')o < 1.4, while the soJid dots with error- 
bars show the bluer clusters in the color range 0.2 < (g'-i')o < 0.8. 
Both sets of points have been background-subtracted as described 
in the previous figure. The solid line is a Gaussian LF with the 
same mean and standard deviation as in the previous figure, match- 
ing the entire GCLF over all colors. 

ones in the opposite sense expected from bolometric cor- 
rections) may be indicators of genuine differences in their 
mass distribution s. In an HST-based survey of several 
nearby ellipticals. lLarsen et al.l {2001) found a mean dif- 
ference AMy = 0.4 mag in the sense that the blue GCLF 
peak was brighter, though large and uncertain differences 
show up from one system to another. They conclude that 
it is unclear whether any true underlying mass differences 
are required to explain this mean difference. 



NGC 3311 has large numbers of both red and blue 
clusters, and using the same numerical technique as de- 
scribed above, we determined their GCLFs separately. 
"Blue" , metal-poor clusters are defined to lie within the 
color range 0.2 < (g' — i')o < 0.8 and "red", metal-rich 
clusters within 0.8 < (g' — i')o < 1.4. The results are 
shown in Figure [5J where we display only the residual 
LFs and not the raw numbers before background sub- 
traction. The blue population falls close to the same 
roughly-Gaussian curve describing the system as a whole. 
For the red population, the small excess of very luminous 
clusters (i' < 23) is evident, but the most striking dif- 
ference is in the range i' = 25 — 26, where a noticeable 
dip is seen just around the nominal turnover point. We 
are, however, reluctant to conclude that this effect is real 
given the uncertainties in the background correction and 
the determination of residuals. A standard Kolmogorov- 
Smirnov two-sample test shows that the two LFs dif- 
fer over the range i < 26 at more than 99% confidence, 
but the strongest statement we feel is justified from the 
present data is that the bright-end slope for the red clus- 
ters is somewhat flatter than for the blue ones. A better 
control field, and perhaps even better selection of GCs 
versus field contamination, will be needed to carry this 
discussion any further. 

4.5. Specific Frequencies 

Finally, we use the radial distributions we have de- 
rived for NGC 3311 and NGC 3309 to estimate the to- 
tal cluster populations around each one, and thus their 
specific frequencies. We use the best-fit Sersic profile 
parameters given above, and integrate them outward to 
large enough radius r max that the GC population to- 
tals converge. For NGC 3309, the particular choice of 
Tmax is quite unimportant, since its profile is steep and 
negligible numbers of GCs are added for r > 3'. How- 
ever, the much more extended profile for NGC 3311 re- 
quires us to integrate outward to r max ~ 9' or about 
140 kpc, beyond which the Sersic model profile adds 
negligible numbers of GCs. This large radius is, how- 
ever, quite consistent with the existing surface p hotom- 
etry (jRichter et al.l 119821 : IVasterberg et al.l 119911 ) which 
shows that NGC 3311 extends detectably to r ~ 150 kpc. 
With these parameters in mind, we find total GC pop- 
ulations brighter than i' = 25.4 equal to A3311 = 7440 
and N 330g = 165. 

Evaluating the actual population totals requires cor- 
recting for the magnitude limit and the fraction of GCs 
fainter than that limit. Assuming a Gaussian-like GCLF 
as discussed in the previous section, with a turnover point 
at My ~ —7.4 ±0.1 and standard deviation ay ~ 1.5, 
then i' = 25.4 is 0.25 ± 0.1 mag or 0.17c brighter than 
the turnover point. This part of the GCLF includes a 
fraction (0.45 ± 0.04) of the total population. We there- 
fore estimate the total GC population for NGC 3311 as 
N t ~ 16500 ± 2000, where the main uncertainty is the 
necessary extrapolation to large radius. 

For NGC 3309 we find N t = 364 ± 210; here, the rel- 
atively large internal error is dominated instead by the 
large uncertainty in the ratio (a^/vei) = 0.080 ± 0.046 
for the two-function fit described above. In this case our 
calculation may be a slight underestimate, not because 
of what may be present at larger radius, but because 
of the uncertain inward extrapolation into the central 
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~ 20". Its GCS core radius is both smaller and consid- 
erably more uncertain (by a factor of two; see above), 
and our profile integration may not have accounted for a 
significant part of the total. The HST photometry to be 
discussed in the next section verifies that more GCs are 
present in its core. 

Their specific frequencies ca n be determined once we 
know their total luminosities. Vas terberg et al.l (|1991l ) 
measure integrated magnitudes Vt = 11.04 for NGC 
3311 and 11.94 for NGC 3309 (we have reconstructed 
these numbers from their Table 5 and their adopted dis- 
tance and reddening). With our adopted (m — M)y = 
33.88, their luminosities are M$ (JV3311) = -22.8 and 
M£(iV3309) = -21.9. The (global) specific frequency 
of the NGC 3311 GCS is then S N = N t ■ io°- 4 ( 15 - M ?) = 
12.5±1.5, and for NGC 3309 SW = 0.63±0.36. It should 
be noted, however, that these total magnitudes are ex- 
trapolations to large radius of the r 1 / 4 profiles that were 
found to fit the inner regions, and for NGC 3311 at least 
a part of the excess cD envelope light would be missed. 
Vasterberg et al. estimate that it should be made ~ 0.4 
mag brighter, which would give M£(A3311) = -23.2. If 
this adjustment is applied, its specific frequency would 
be lowered to Sn = 8.7±1.0. However, it is not clear that 
the luminosity should be increased without also making a 
similar adjustment to the total cluster population, since 
we might have missed some of the GCs by using only 
the extrapolation from the inner profile. For the present 
time, and pending a more quantitative measurement of 
the cD envelope light, we use Sn — 12, while Sn ~ 9 
should be considered an extreme lower limit. 

By contrast, NGC 3309 has a remarkably low deduced 
cluster population, even with the uncertainties about the 
core population. Giant E galaxies in rich cluster environ- 
ments wi th Sn 5, 1 are rare, but not unprecedented. For 
example, iBaum et al.l (|1995j ) found Sn < 1 for one of 
the Coma cluster giants NGC 4881. It is tempting to 
imagine that an object like NGC 3309 might have had 
much of its original GCS removed by tidal stripping in 
passages through the cluster core. Regardless of its his- 
tory, it is clear that today, NGC 3311 completely dom- 
inates the GCS population in the core Hydra field, in 
equal measure because of its higher total luminosity and 
its higher intrinsic Sn- The wide ranges of Sn that we 
find among these giant and structurally simple E galaxies 
are still not fully understood and may require combina- 
tions of formation conditions, later dynamical evolution 
in the potential field of their cluster, and a ccounting for 
the presence or ab sence of hot halo gas (e.g. lHarrial2001l : 
lMcLaughlinlfT999l) . 

As a final test of the relative GC populations of the 
two galaxies, we ran another solution for Eq. [4] where we 
imposed the assumption that the total number of GCs 
in each galaxy was in proportion to the luminosity of the 
galaxy; that is, they have the same Sn- This requirement 
converts to A^ G c(3311)/iV G c(3309) ~ 2.5. The best-fit 
solution gave x^(mirt) ~ 2.8, worse than our optimum 
solutions described above (xt = 2.0) though not dramat- 
ically so. The implication is simply that the parameters 
for the NGC 3309 system are poorly constrained within 
broad limits, a statement that is already reflected in the 
large uncertainty on its deduced best-fit cluster popula- 
tion. 
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Fig. 9.— WFPC2 fields in the Hydra core. The panel on the left 
shows the data for the field in which the PCI chip is centered on 
NGC 3311, while the right panel shows the field centered on NGC 
3309. All measured objects brighter than I = 25.0 are plotted. 
The xy positions on the frames are plotted in arcseconds. Notice 
that the orientations of the two fields are different: for NGC 3311, 
North is 109.8° counterclockwise from vertical, and for NGC 3309, 
North is 19.7° counterclockwise from vertical. Within the "NGC 
3309" field, the center of NGC 3311 lies to the lower left; notice 
the increase in numbers of measured points there. 

NGC 3311 clearly holds one of the very largest globular 
cluster systems in the local universe, rivalling that of the 
classic high— S at prototype M87 which has N t ~ 13000 
and S N ~ 13 ([McLaughlin et all 11994 iTamura et ail 
12006ft . The NGC 3311 GCS is spatially very extended, 
marked by a huge core radius c = 94" = 24 kpc or effec- 
tive radius r e — 177 " = 45 kpc also co mparable to that 
of the GCS in M87 (|McLaughlin|[T999T) . 

5. WFPC2 PHOTOMETRY 

To gain additional insight into the relative properties 
of the GCSs around NGC 3309 and 3311, we employed 
imaging in (V, I) from the HST WFPC2 camera that was 
available in the HST Archive. Two WFPC2 fields were 
taken as part of GO program 6554 (PI Brodie) in Cycle 6, 
each with total exposures of 3700 sec in F55W and 3800 
sec in F81AW. In one field, the PCI chip was centered 
on NGC 3311, and in the second field, the PCI chip was 
centered on NGC 3309. Photo metry of the NGC 3311- 
centered field was published in iBrodie et al.l (|2000f ) , but 
to our knowledge, no measurements from the NGC 3309 
imaging have been published. 

We have carried out photometry on both of these 
WFPC2 fields with daophot in the same manner as de- 
scribed above. A typical globular cluster with a half-light 
radius of ~ 3 pc will, at the distance of Hydra, have an 
intrinsic FWHM < 0"03, making it entirely feasible to 
perform PSF-fitting photometry with the WFPC2 reso- 
lution of 0"1. For each of the'four WFPC2 CCDs, we 
registered and co-added the individual exposures down- 
loaded from the HST Archive, and performed the usual 
sequence of daophot FIND/PHOT/PSF/ALLSTAR. De- 
tected objects on the V and / frames that matched loca- 
tions within 2 pixels, and with ALLSTAR goodness-of-fit 
Xv.i < 2.0, were kept for the final data l ist. Calibrations 
follow ed the normal prescriptions in Holtzman et al. 
(1995), based on 0V5 aperture photometry of the bright- 
est objects. The data for the four CCDs were then com- 
bined to yield the results shown in Figures l9l and [T0l The 
photometric limit near I ~ 25 or i' ~ 25.3 is about half 
a magnitude shallower than for our GMOS data. 

Our PSF-based photometry for NGC 3311 has gen- 
erated a sl i ghtly tighter CMD than the work of 
IBrodie et al.l (|2000f ) (who used only aperture photome- 
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Fig. 10. — Photometry for the measured starlike objects in the 
two WFPC2 fields shown in the previous figure. 
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Fig. 11. — Photometry for measured objects within 30" of the 
center of NGC 3309. 

try), but no major differences. Though the (V — I) color 
index is notoriously insensitive to metallicity (with an 
intrinsic difference of only 0.2 mag between the red and 
blue seq uences), a clear trac e of bimodality can be seen 
(also see iBrodie et al.l [2000) . Notably, on the red side 
at (V — I) ~ 1.2, the upward extension of the GC se- 
quence into the high-mass UCD-like range that is more 
evident in our more com prehensive GMOS photometry 
(jWehner fc Harris! [200l can also be seen here. In gen- 
eral, both CMDs in Fig. [TOl resemble each other closely. 

In the xy plots of Fig. El we see the "NGC 3309 prob- 
lem" once again. In the so-called NGC 3309 field, a pop- 
ulation of GCs is present on and near the PCI chip cen- 
tered on 3309; but the presence of the NGC 3311 GCS 
is extensive everywhere over the field (note particularly 
the lower left region of the xy plot where the increase in 
GC number density toward the center of NGC 3311 is 
particularly obvious). Even though the exposures were 
intended for study of the NGC 3309 system, the ma- 
jority of detected objects in it actually belong to NGC 
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Fig. 12.— Color-magnitude data for both WFPC2 fields com- 
bined, excluding any objects within 40" of the center of NGC 3309 
(within the dashed circle in Fig.|9]l. This CMD is completely dom- 
inated by GCs belonging to NGC 3311. 

3311. From our quantitative double-profile solutions in 
the previous section, it can easily be shown that for any 
location further than 40" from the center of NGC 3309, 
the GCS of NGC 3309 contributes less than 10% of the 
GC counts. This transition boundary beyond which the 
NGC 3309 system becomes unimportant is indicated by 
the dashed circle in Fig. [5] The CMD for this inner zone 
is shown in Figure [TT1 where we have included the ~ 360 
detected objects within an even smaller radius of 30". 
In this diagram, which gives us our best look at the GC 
population of NGC 3309 alone, the high-luminosity end 
of the CMD is not as well populated as in the NGC 3311 
GCS (perhaps a simple result of small-number statistics), 
but otherwise there are no striking differences to note. 

In Figure [T2l we have combined the photometry from 
both WFPC2 fields, excluding only the points within the 
40" zone around NGC 3309. This CMD, dominated by 
NGC 3311, will be analyzed below simultaneously with 
the GMOS data, for bimodality characteristics of the 
metallicity distribution. 

6. THE METALLICITY DISTRIBUTION AND ANALYSIS OF 
THE BIMODALITY 

A primary goal of our study was to investigate the mul- 
timodal structure of the GC color and metallicity distri- 
bution, and in particular to find out whether or not the 
blue GC sequence would reveal the same mass-metallicity 
relation as in other large galaxies. Having measured sev- 
eral other properties of the system, including the fact 
that the GC population in our target field is dominated 
by NGC 3311, we now discuss the color distribution. 

The histogram of colors (g' — i') over any selected mag- 
nitude range drawn from Fig. [2] indicates immediately 
that the color distribution is broad and not well fit, for 
example, by a sing le Gaussian-type curve (also found by 
IBrodie et all 120001 from their WFPC2 photometry). A 
better description of the system is by conventional "blue" 
and "red" sequences whose color distributions partially 
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overlap. The bright-end extension of the red sequence 
in part icular reaches upward i nto the UCD luminosity 
range (|Wehner fc Harris! I2007T ) . We therefore experi- 
mented with bimodal fits in a variety of ways, to help 
define the mean colors and dispersions of these two com- 
ponents. 

To carry out quantitative fits to the color distribution, 
we e mployed the statistical code RMIX (|MacDonaldl 
2007), a library written with the statistics programming 
language, R. This package allows the user to fit data with 
various functional shapes and to set constraints on them 
at a range of possible levels. RMIX does not restrict 
the fitting to unimodal or bimodal forms; the number of 
components adopted to fit the data is user-defined, as are 
the functional forms of the curve (e.g. Gaussian, Poisson, 
binomial) and the starting values for the modes and dis- 
persions of each component (e.g. the peak and standard 
deviation, in the Gaussian case). In addition, the disper- 
sions can be constrained to be equal (homoscedastic) or 
allowed to differ between components (heteroscedastic), 
as appropriate. The numerical values of the dispersions 
can also be set to previously determined numbers, or can 
be determined by the fit, in any combination between the 
components. 4 

The KMM soft ware (jMcLachlan fc Basfordl 119871 : 
lAshman et alJ[l994h is a well known code frequently used 
in the literature for bimodal fits to GCS color distri- 
butions (see, e.g. lLarsen et al.l 120011 ; iMieske et al.ll2006l 
among many such examples). We have experimented 
with both KMM and RMIX in this study. Both give 
quite similar results for a given dataset, but we find that 
the RMIX package offers more flexibility in its choice of 
fitting functions as well as better flexibility in managing 
user-defined constraints on the fitting parameters. It is 
also more robust, successfully converging on multimodal 
fits in small-sample cases where KMM would not run. 
We recommend RMIX as an attractive and statistically 
rigorous alternative. 

Figure [13] shows some of our results for four different 
magnitude subsets of our (g',i') data. For each mag- 
nitude bin shown, we fit a pair of Gaussian curves to 
the (g' — i') histogram. Reasonable initial guesses for 
the peaks and dispersions were provided, and the com- 
ponent means and dispersions were determined by the 
fit itself with the condition that the dispersions of the 
two Gaussians were required to be the same. The re- 
sults are summarized in Table [2] Here, the first column 
gives the magnitude range; column 2 the deduced pro- 
portions of the blue and red components; columns 3 and 

4 the mean colors of the blue and red sequences; column 

5 the standard deviation of each sequence; and the last 
column the goodness of fit. Bimodal fits are strongly 
preferred over unimodal ones for the two brighter bins. 
In the two faintest bins the overlap between the red and 
blue components increases because of the combination of 
field contamination and larger photometric uncertainty, 
and the distinction between bimodal and unimodal fits 
becomes less clear. 

In Fig. 021 one can see that the peak of the 

4 The complete code, available for a variety of plat- 
forms, is publicly available from Peter MacDonald's website at 
http://www.math.mcmaster.ca/peter/mix/mix.html The same 
site gives links to an extensive bibliography with turther examples 
of its use. 
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Fig. 13. — Color distribution for the globular cluster popula- 
tion in NGC 3311, as a function of magnitude. Number of GCs 
per unit (</ — i')o interval is plotted for four different luminosity 
bins as listed in Table 2. The interval i' = 22 — 23 is at upper 
left, 23 - 24 at upper right, 24 - 25 at lower left, and 25 - 26 at 
lower right. In each panel, the bimodal Gaussian fit as determined 
by the RMIX code is shown overplotted on the histogram of the 
data. The two individual Gaussians matching the metal-poor and 
metal-rich components are shown by the red lines, while the green 
envelope is the direct sum of the two components. Small filled tri- 
angles on the x-axis of each plot mark the mean color (Gaussian 
peak position) of each component. Note that the metal-poor com- 
ponent becomes progressively redder at higher luminosity, while no 
significant change occurs for the metal-rich component. 

bluer component shifts closer to the red component 
at brighter magnitudes, suggesting the presence of a 
mass-metallicity trend for at least one of the modes. 
To explore this further, we separated the GC popula- 
tion into blue and red halves by taking the intervals 
0.3 < (g' - i') < 0.8 (blue) and 0.8 < {g 1 - i') < 1.2 
(red), and then calculating the mean color of each mode 
in half-magnitude bins. The combined results can be 
seen in Figure 1141 Linear least-squares fits over the en- 
tire GC luminosity range 22 < i' < 26, weighted by 
1/er 2 where the uncertainty of each mean color is a, give 
A(g'-i')/Ai' = 0.044±0.011 for the blue sequence with 
a correlation coefficient of the fit equal to 0.995, indicat- 
ing that a MMR is present at high significance. For the 
red sequence, A(g' - i')/Ai' = 0.012 ± 0.021, indicating 
no significant trend with magnitude. In Fig. [14] we show 
it simply as an average (<?' — i')o,red = 0.94. 

Field contamination at magnitude levels i' < 25 is 
unimportant, so at these relatively bright levels in the 
color-magnitude diagram the existence of an MMR seems 
secure. Fainter than this, however, contamination be- 
comes progressively worse and is likely to affect these 
simple linear fits to uncertain degrees. Thus for the 
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Fig. 14. — Color-magnitude diagram for the NGC 3311 globular 
cluster system, plotted as absolute magnitude versus dered- 
dened color. Mean points for the blue and red globular cluster 
sequences are overplotted, derived as described in the text. The 
error bars on each point are the one-sigma uncertainties in (<?'— i')o- 
The best-fit mean lines determined by linear least squares are 
drawn through each sequence. For the blue sequence, the slope is 
A(cj' — = 0.044 ± 0.011, corresponding to a heavy-element 

abundance scaling with cluster mass of Z ~ M 0,6 . The red GC 
sequence shows no net change in color with magnitude, and is as- 
sumed for the purposes of this plot (see text and previous figure) 
to be at a fixed color of (g' — i')o = 0.94. 

fainter parts of the sequences (roughly, > —9.5) 
we cannot test with any confidence whether or not the 
blue sequence really does continu e to get "bluer" as it 
goes down. In lHarris et al.l (|2006f ). we claimed that the 
blue sequence became approximately vertical (that is, 
no MMR) for luminosity levels Mj > -9.5. That claim 
was based on bimodal fits to the color distributions in 
finely spaced magnitude bins for the eight BCGs stud- 
ied there; however, in most of those galaxies the issue 
of faint-end field contaminatio n is also a concern just 
as it is here. For NGC 4594, ISpitler et all (|2006f) sug- 
gest that the blue-sequence MMR continues downward 
roughly linearly to fainter luminosities, though their two 
faintest blue-sequence bins are not distinguishably dif- 
ferent in mean color, and their conclusion of the same 
MMR at all luminosities is in part forced by their as- 
sumption of a linear fit. To quote them directly, "Indeed, 
whether the faintest bins follow the trends established by 
the brighter points is indeterminate." For our NGC 3311 
data, removing the two faintest bins, where we expect the 
contamination to be most significant, as well as the top 
bin, which m ay be affected by the pre sence of Ultra Com- 
pact Dwarfs (|Wehner fc Harrisll2007D , yields a similar fit 
of A{g' -i')/Ai' = 0.040 ±0.011. As such, the resulting 
slope does not appear to be driven by contamina t ion in 
the faintest magnitude bins. Lastly, Mics ke et al.l (2006) 
use the combined GCS data for 79 Virgo galaxies to test 
for color trends by combining their target galaxies into 
four groups by luminosity. In their data, field contami- 
nation is much less of a concern because individual GCs 
were carefully selected by a series of morphological and 
photometric criteria. For the highest-luminosity galaxies 
(M87, M49, and M60; see their Figures 1 and 2) it can 



reasonably be claimed that the blue-sequence MMR con- 
tinues roughly linearly over the entire luminosity range 
— 11 > M z i > —7. For the other groups containing a 
total of 76 galaxies from giants to dwarfs, the blue GC 
sequence clearly shows a color trend at bright levels but 
could reasonably be claimed to become nearly vertical 
for M z i > —9, corresponding roughly to My > —8.5. 

At the present stage in this subject, we cannot yet say 
on the basis of the empirical evidence whether or not 
the blue-sequence MMR persists at all luminosity levels. 
More complete data allowing for accurate removal of field 
contamination will be needed to trace it out at the crucial 
faint levels. All studies to date agree, however, that a 
MMR exists for the blue clusters brighter than about 
Mi ~ —9.5, and that this seems to be a nearly universal 
phenomenon. 

Comparing our MMR for NGC 3311 with those found 
previously for other large galaxies requires transfor- 
mations between photometric system s, si nce H06 used 
B - I), whereas IMieske et al.l (|2006f ) and lStrader et all 
20061) used (g' - z'). Ideally, we would like to convert 
(g 1 — i') to cither one of these other indices, but transfor- 
mations based directly on observations of standard glob- 
ular clusters (the Milky Way clusters, for example) do 
not yet exist. For the purposes of a preliminary compar- 
ison, we choose instead to transform the observed blue- 
sequence slope A(g' — i') into a metallicity scale A[Fe/H] 
through the use of models. Representative calc ulations of 
integr ated c olors for SSP systems are given bv lMarastonl 
(2005) and iGirardi et all (|2004f ). For the range with 
[m/H] < —0.5, these two sets of models agree well with 
each other and give A{g' - i')/A[m/H] = 0.18 ± 0.03 for 
the metallicity regime we are concerned with, i.e. [m/H] 
< —1.2. Then, converting color index to metallicity, and 
magnitude to luminosity and hence mass, we find that 
the MMR for the blue sequence in NGC 3311 corresponds 
to a heavy-element abundance scaling with clust e r mas s 
Z ~ M°- 6±0 - 2 . As we did in IWehner fc Harris! (f2007h . 
we assume that (M/L) does not change systematically 
with luminosity in order to convert L to M and thus 
derive the scaling. By comparison, H06 found a mean 
scal ing Z ~ M 0,5 5 for t heir sample of eight giant galax- 
ies. IMieske et all £2006) found Z ~ M°- 48±0 08 for their 
Virgo galaxy sample. Our results for NGC 3311 are en- 
tirely consistent with the other studies within measure- 
ment uncertainties. 

As an additional check of the (g 1 — i 1 ) analysis, we 
perform the same kinds of bimodal Gaussian fits to the 
WFPC2 data in (V-I). Here, the smaller size of the (V- 
I) sample, and its lower sensitivity to metallicity, prevent 
us from tracing out mean color trends with magnitude 
with any confidence, so we use a color histogram defined 
from only a single magnitude range 21.5 < I < 24.5 over 
which the field contamination and photometric scatter 
are minimized. As before, we exclude the region within 
40" of NGC 3309. 

The histogram of 1730 such objects is shown in Fig- 
ure [HI We ran RMIX under similar assumptions as de- 
scribed above (two Gaussians, with means, dispersions, 
and proportions determined by the fit). For the blue, 
metal-poor component we find (V — I) = 1.10 ± 0.07, 
ay-i = 0.12 ± 0.02, and proportion of the total pop- 
ulation piblue) = 0.58 ± 0.39. For the red, metal-rich 
component, (V - I) = 1.28 ± 0.07, <7y_/ = 0.11 ± 0.02, 
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Fig. 15. — Color distribution for globular clusters in NGC 3311, 
from the WFPC2 (V — I) photometry discussed in the text. No 
dereddening has been applied. The double-Gaussian bimodal fit to 
the data is shown, where the solid curve equals the sum of the two 
individual components shown in the dashed lines. The parameters 
for the two components are listed in the text. 

and p(red) = 0.42 ± 0.39, with an overall \v = 1-22. 
The relative proportions of the two components are quite 
uncertain - again, a direct result of the insensitivity of 
V — I to metallicity and thus the heavy overlap between 
the two components - but agree nominally with the near- 
equal proportions we found from the (<?' — i') solutions 
(Fig.rr3J). For comparison purposes, we also tried a fit of 
a single Gaussian curve to the entire data: this yielded a 
best-fitting mean (V — I) = 1.175 ± 0.004 and dispersion 
a V -i = 0.146 ± 0.003, with goodness-of-fit X v = 1-47. 
The bimodal fit performs better at matching the data. 

Although (V — I) is less sensitive to metallicity, it has 
the advantage over the SDSS indices that its conversion 
to [Fe/H] is more well determined. With Ey-i = 0.11, 
we have (V-I) (blue) = 0.99±0.07and (V-I) (red) = 
1.17±0 .07. Adopting a conversion (V-I)p = .16 [Fe/H] 
+ 1.15 (|Harris et al.ll2000HBarmbv et alJ 12000) based on 
calibration from the Milky Way and M31 clusters then 
gives [Fe/H] ~ —1.2 for the blue sequence and ~ +0.2 
for the red sequence. These mean metallicities should 
be taken only as rough estimates, since the red-sequence 
mode in particular requires an extension of the conver- 
sion equations into a high-metallicity regime that the 
calibrating GCs in M31 and the Milky Way do not ad- 
equately cover. However, the mean (V — I)q values 
themselves are close to what would be predicted from 
the established correlations with host g alaxy luminosity 
(|Larsen et al.ll2001t iStrader et aTll2004ti . namely 

(V-I)(blue) = -0.011M£ + 0.71, (10) 
(V-I)(red) = -0.018My + 0.80 (11) 

For (iVGC3311) = -22.8 these relations predict 
(V — I)o = 0.96 and 1.19 for the blue and red modes 
respectively, in good agreement with our observed mean 
colors. Supergiant ellipticals such as this one are the 
most logical places to search for GCs of super-Solar 
metallicity, although higher quality spectra than cur- 
rently exist will be required to establish their abundances 
accurately. 



7. SUMMARY 

We have presented an extensive new photometric study 
of the large globular cluster system around NGC 3311, 
the central cD galaxy in the Hydra cluster. The relatively 
wide field and deep limiting magnitude of our data from 
the Gemini/GMOS camera has allowed us to analyze the 
GCS with a sample of several thousand globular clus- 
ters; our data reach to i'(lim) ~ 26.3, just beyond the 
turnover point of the GC luminosity function. Among 
our findings are these: 

• The classic metal-poor blue GC sequence and the 
metal-richer red sequence are both present at their 
expected colors, and with nearly equal populations. 
We show the color distributions both from our 
GMOS (g'-i') photometry and from HST WFPC2 
archival data in (V — I). 

• Bimodal fits to the complete color distribution, 
subdivided by magnitude, confirm that the blue 
sequence has the same correlation of progressively 
increasing metallicity with GC mass that was pre- 
viously found in many other large galaxies; the 
correlation we find corresponds to a scaling of GC 
metallicity with mass of Z ~ M 6 . 

• By contrast, the red sequence shows no change of 
mean metallicity with mass, though it does extend 
upward to much higher than normal luminosity 
into the UCD-like range, strengthening the poten- 
tial connections between massive GCs and UCDs. 

• Confirming previous but much more preliminary 
investigations, we find that the GCS in NGC 3311 
falls definitively into the "high specific frequency" 
category inhabited almost uniquely by supergiant 
ellipticals (cD's or BCG's) at the centers of rich 
clusters. NGC 3311 has an estimated total popu- 
lation of about 16,000 clusters and a specific fre- 
quency Sn — 12. 

• At the core of the Hydra cluster, another giant 
elliptical NGC 3309 is sitting projected just 100" 
from the cD NGC 3311. We use our database to 
solve simultaneously for the spatial structure and 
total GC populations of both galaxies at once. The 
results show that NGC 3309 contributes only a few 
percent of the GC population over the entire field, 
and has an unusually (though not unprecedentedly) 
low value 5^(^3309) = 0.6+0.4. The central cD is 
completely dominant in the central Hydra region. 

• The GC luminosity function, which we measure 
down to the "turnover" point at Mj ~ —8.4, also 
has a normal structure relative to other giant ellip- 
ticals. 

Additional observations that could make clear im- 
provements to the understanding of the NGC 3311 sys- 
tem would include imaging with a still wider field (to 
establish a true "control" field for background measure- 
ment), and multicolor photometry with still more sensi- 
tive color indices to improve the definition of the metal- 
licity distribution function. The resources held by this 
enormous and fascinating GC system are far from being 
exhausted. 
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